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We present the first measurement of the Q 2 -dependence of the neutron spin structure function 
52 at five kinematic points covering 0.57 (GeV/c) 2 < Q 2 < 1.34 (GeV/c) 2 at x ~ 0.2. Though the 
naive quark-parton model predicts <?2 = 0, non-zero values occur in more realistic models of the 
nucleon which include quark-gluon correlations, finite quark masses or orbital angular momentum. 
When scattering from a non-interacting quark, can be predicted using next-to-leading order fits 
to world data for g™. Deviations from this prediction provide an opportunity to examine QCD 
dynamics in nucleon structure. Our results show a positive deviation from this prediction at lower 
Q 2 , indicating that contributions such as quark-gluon interactions may be important. Precision 
data obtained for g" are consistent with next-to-leading order fits to world data. 

PACS numbers: 13.60.Hb, 14.20.Dh, 13.88.+e, 24.85. +p 



Over the past 30 years, significant progress has been 
made in understanding the spin structure of the nucleon 
through measurements using polarized deep- inelastic lep- 
ton scattering (DIS). Most of these experiments were 
focused on precise measurements of the spin structure 
function g\ . In the naive quark-parton model (QPM) , g\ 
is directly related to the contributions of the individual 
quark flavors to the overall spin of the nucleon (see e.g. 
Ref. 0]). Sum rules based on this simple model have 
provided fertile ground for understanding the origin of 
the nucleon spin in terms of quark degrees of freedom. 
In addition, next-to-leading-order (NLO) analyses of the 
world gi data (see e.g. Refs. 0)0) have provided indi- 



rect information about the role of gluons in the nucleon's 
spin. 

The QPM is expected to be valid in the scaling limit, 
where the four-momentum transfer squared — Q 2 and en- 
ergy transferred v approach infinity. As Q 2 becomes 
large, the electron-nucleon interaction can be described 
by scattering from a massless, non-interacting quark car- 
rying a fraction x = Q 2 /2Mv of the nucleon's momen- 
tum (M is the nucleon mass). At finite Q 2 and v, ef- 
fects such as gluon Bremsstrahlung, vacuum polarization 
and vertex corrections can be accurately calculated us- 
ing perturbative QCD. For g\, non-perturbative QCD 
processes, such as quark-gluon and quark-quark correla- 
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tions, are suppressed relative to the asymptotically-free 
contributions by factors of 1/Q and l/Q 2 , respectively. 

Polarized DIS also provides information about a sec- 
ond spin structure function, g 2 , which is identically zero 
in the QPM |4|. Interest in g 2 arises because, unlike 
gi , contributions from certain non-perturbative processes 
enter at the same order in Q 2 as the asymptotically-free 
contributions. An appropriate formalism for understand- 
ing g 2 is the operator product expansion (OPE) || Q 
which is a model-independent approach based directly 
on QCD. Here, the unknown hadronic currents relevant 
for polarized DIS are expanded in terms of quark and 
gluon operators and grouped by factors of (1/Q) T ~ 2 , 
where r = 2,3,4,... is known as the twist of the op- 
erator. Both c/i and g 2 contain a contribution from a 
twist-2 operator that corresponds to scattering from a 
massless, non-interacting quark. Operators with higher- 
twist represent contributions from non-perturbative pro- 
cesses such as quark-quark and quark-gluon correlations, 
and from quark mass effects. Because these correlations 
are responsible for quark confinement, higher- twist ef- 
fects must be included in any realistic model of the nu- 
cleon 0. 

Using the fact that g\ and g 2 contain the same twist-2 
operator, Wandzura and Wilczek p| derived the following 
expression for the asymptotically- free contribution to g 2 , 
in terms of g±, 



g 2 vw (x,Q 2 ) = -g 1 (x,Q 2 ) + f Q ' } dx. 



(1) 



The world data for g\ cover a broad range in x and Q 2 
with relatively high precision. Polarized parton distribu- 
tions can be extracted from NLO fits to the data and 
evolved to different values of Q 2 using the Dokshitzer- 
Gribov-Lipatov-Alarelli-Paresi (DGLAP) procedure 0, 
El EH EH • Because the world data are at relatively large 
Q 2 , where higher- twist effects should be negligible, these 
evolved parton distributions allow one to calculate the 
twist-2 part of g\, and therefore g 2 vw , in most kinematic 
regions accessible today. Precise measurements of g 2 at 
specific values of x and Q 2 can be compared to g 2 , 
providing a unique opportunity to cleanly isolate higher- 
twist contributions. 

Previous measurements of g 2 [l2L Hij , HiL Il5| were 
aimed at testing OPE sum rules and covered a wide 
range in x, at an average Q 2 of 5 (GeV/c) 2 . Proton 
data show general agreement with g 2 rw but lack the pre- 
cision needed to make a definitive statement about the 
size of higher-twist effects by direct comparison. Neu- 
tron results have much larger uncertainties and cannot 
distinguish between g\ vw and g 2 — 0. We present a new 
measurement of the Q 2 -dependence of 172 for the neu- 
tron at low Q 2 , while keeping x approximately constant. 
With statistical uncertainties more than an order of mag- 
nitude smaller than existing data, these results allow, for 



the first time, a precise comparison with g 2 vw to study 
the Q 2 -dependence of higher-twist effects. 

Longitudinally-polarized electrons were scattered from 
a polarized 3 He target in Hall A 0] at the Thomas 
Jefferson National Accelerator Facility (Jefferson Lab) 

through the inclusive process 3 He(e, e'). We measured 
g 2 at five values of Q 2 between 0.57 (GeV/c) 2 and 1.34 
(GeV/c) 2 , with x ~ 0.2. The kinematics for this x are 
well-matched to the beam energies available at Jeffer- 
son Lab and are in a region where g^ w is relatively 
large. Low Q 2 provides access to the region where higher- 
twist effects are expected to become important. The 
invariant mass squared of the photon-nucleon system, 
W 2 = M 2 + 2Mv - Q 2 , was kept > 3.8 (GeV/c) 2 to 
minimize resonance contributions. The spin of the 3 He 
nuclei could be oriented either longitudinal or transverse 
to the incoming electron spin. 

The parallel and perpendicular cross-section differ- 
ences, Act 1 1 and Actj_, are defined in terms of helicity- 
dependent cross-sections o % * as, 



Actii 



Aft 



Act 



(2) 



where the single and double arrows refer to the beam 
and target spin directions, respectively, relative to the 
incident electron momentum direction. An up (down) 
arrow refers to spin aligned parallel (anti-parallel) to the 
momentum. The double right arrow refers to target po- 
larization in the scattering plane, perpendicular to the 
incident electron momentum, with the scattered electron 
detected on the side of the beam line towards which the 
target spin points. These cross-section differences are re- 
lated to g\ and g 2 by [lj, 



Actm 



Act 1 = 



Aa 2 E' 
Q 2 EM v 
2xMg 2 ( Xl Q 2 )] 



[{E + E' cose)g 1 (x,Q 2 ) 



4,a 2 E' 
Q 2 EM v 



E' sme[ gi (x,Q 2 ) 



AxEM 



g 2 (x,Q 2 )], 



(3) 



(4) 



where E and E' are the incident and scattered electron 
energies and 9 is the laboratory scattering angle. 

Polarized electrons were produced by photo-emission 
from a strained GaAs crystal using circularly polarized 
laser light. The average beam polarization was Pb = 
(76 ± 3)% as measured using both M0ller and Compton 
polarimeters |l6| and the average beam current was 12.0 
/j A, with an uncertainty of ±1%. The helicity of the 
electron beam was flipped on a pseudo-random basis at a 
rate of 30 Hz to minimize helicity-correlated systematic 
effects. The beam energy ranged from 3.5 GeV to 5.7 
GeV and was measured with a relative uncertainty below 
10~ 3 using both elastic electron-proton scattering and a 
magnetic field measurement. 
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Scattered particles were detected by either one of a 
pair of magnetic spectrometers arranged symmetrically 
on either side of the beam line 0] . The momentum of 
the particles was determined by reconstructing their tra- 
jectories using drift chambers. Electrons were identified 
using a threshold gas Cherenkov detector and a two-layer 
lead-glass electromagnetic calorimeter. The largest back- 
ground was from ir~ photo-production in the target. The 
ratio of pion rate to electron rate detected in the spec- 
trometer was < 3.3. The ratio of the pion asymmetry to 
the electron asymmetry was < 4.1. The pion rejection 
factor was > 10 4 : 1 with an electron detection efficiency 
of > 97%, which was sufficient to reduce the pion con- 
tamination to a negligible level. Electrons scattered in 
the entrance and exit windows of the target cell were 
removed using software cuts. 

To study polarized neutrons, a target containing 3 He 
nuclei was polarized using spin-exchange optical pump- 
ing 01- The 3 He ground state is dominated by the S- 
state in which the two proton spins are anti-aligned and 
the spin of the nucleus is carried entirely by the neu- 
tron. The 3 He is contained in a sealed, two-chambered, 
aluminosilicate glass cell, along with a small quantity of 
N2 and Rb to aid in the polarization process. Polar- 
ized 3 He is produced in the spherical upper chamber by 
first polarizing Rb atoms with optical pumping. These 
atoms can transfer their spin to the 3 He nucleus during 
binary collisions. Incident electrons scatter from the po- 
larized 3 He in the cylindrical lower chamber which is 40 
cm long with side walls of thickness » 1.0 mm and end 
windows of thickness 130 /im. The 3 He density as seen 
by the beam is 2.9 x 10 20 /cm 3 . The average in-beam 
target polarization was P t = (40.0 ± 1.4)% as measured 
using both nuclear magnetic resonance |l8( and electron 
paramagnetic resonance |l9| . 

The 3 He spin structure functions g\ and 172 were ob- 
tained using our measured cross-section differences and 
Eqs. J3J and (@J. The ratios of these cross-section dif- 
ferences to the total unpolarized cross-section are de- 
fined as the longitudinal and transverse asymmetries, 
respectively. To achieve our desired statistical preci- 
sion on <?2 1 the raw (uncorrected) transverse asymmetries 
were measured to a statistical uncertainty of ~ 10~ 4 . 
It was also necessary to keep sources of false asymme- 
tries at, or below, this level. A feedback system was 
used to keep helicity-dependent beam charge asymmetry 
below 50 x 10 -6 for a typical run. Quasi-elastic polar- 
ized electron scattering from an unpolarized 12 C target 
was used to measure the false asymmetry and yielded 
(67 ± 46) x 10~ 6 . Misalignment of the target polariza- 
tion direction was typically less than ±0.3°, implying a 
negligible contribution from Aery when measuring Aa±. 
For each kinematic setting, equal quantities of data were 
collected with the overall sign of the electron helicitics 
reversed and/or the target polarization rotated by 180°. 
The absolute value of the measured asymmetries for each 



of the beam helicity and target spin combinations were 
consistent with each other within the statistical uncer- 
tainties, indicating that there were no significant false 
asymmetries. 

The raw cross-section differences were determined with 
a relative uncertainty of 6% and were checked against 
clastic data taken during this experiment. Spectrome- 
ter acceptances were provided by the Hall A Single- Arm 
Monte Carlo simulatio n polEU used for previous polar- 
ized 3 He experiments |22l |23| . Corrections were applied 
for beam and target polarizations, P& and P t , and for ra- 
diative effects using a modified version of the code POL- 
RAD 2.0 E-pl for internal corrections, and the formalism 
in Refs. [25L |26J for external corrections. Input to the 
radiative corrections came from fits to measured polar- 
ized cross-sections in the quasi-elastic region [27j, fits to 
measured and q" d ata in the resonance region j2^ 
and a fit to 5™ /F™ [23j, along with calculations of g^ /w , 
in the deep-inelastic region. Our measured data were 
used to guide the fits near x ~ 0.2, and the systematic 
uncertainty in the radiative corrections is dominated by 
uncertainties in the fits where data are sparse. 

To obtain neutron spin structure functions, a correc- 
tion was applied to the measured values of g 1 
based on a model 3 He wavefunction 128 



3Ho and ff 2 Hc 



1 



91,2 



where P„ 



P n + 0.056 



3 He 
01,2 



(0.014 -2P p )gl, 



-0.028l°'° 3 ^ and P n 



0.86 



+0.009 
-0.004 



(5) 



are the 



effective proton and neutron polarizations in 3 He |28tl2£ 
l30| . For calculations of g\ we used the average of scenar- 
ios 1 and 2 in the Blumlein and Bottcher (BB) NLO fit 
to world data 0, evolved to our Q 2 using the DGLAP 



procedure d, l9ulOL Hlf . For g\ , values for g\ 
calculated using BB and Eq. 



WW 



were 



TABLE I: Results for an d <?2 from this experiment. The 
two uncertainties are statistical and systematic, respectively. 



Q' 2 


X 


E 


9i 


92 


(GeV/c) 2 




(GeV) 


(1(T 3 ) 


(10- 3 ) 


0.57 


0.16 


3.465 


-71.9 ±5.9 ±8.3 


78.3 ± 9.0 ± 11.0 


0.78 


0.18 


4.598 


-55.7 ±7.5 ±7.2 


75.2 ±13.0 ±10.1 


0.94 


0.19 


4.598 


-51.2 ± 5.9 ± 7.4 


67.5 ±9.0 ±7.2 


1.13 


0.19 


5.727 


-42.6 ±5.3 ±7.0 


48.3 ±8.9 ±6.0 


1.34 


0.20 


5.727 


-36.9 ±5.5 ±6.8 


55.0 ±9.7 ±6.8 



Results for g™ and g% from this experiment are given 
in Table [I] The largest contributions to the systematic 
uncertainties in and g% come from the uncertainties in 
the cross-section differences, and, for g^ , from uncertain- 
ties in the models for the transversely polarized quasi- 
elastic and resonance-region tails used for radiative cor- 
rections. 

Figure ^ shows our data for as a function of Q 2 . 
Also shown are g™ predictions from the NLO analyses of 
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BB and M. Hairi et a/.(AAC03e) Q. The uncertainty in 
the BB curve comes from propagating the uncertainties 
in the parton distribution functions used in the fit. The 
world data for at our x are at significantly larger Q 2 
where higher- twist contributions should be suppressed. 
The good agreement with these evolved fits indicates 
that, within the uncertainties, we are not seeing higher- 
twist effects in g™. This also gives confidence that we 
can use the BB fit to calculate gY W without introducing 
significant higher-twist effects from g\. 
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FIG. 1: Results for g" as a function of Q 2 . Data are shown 
with statistical uncertainties only, with the systematic un- 
certainty indicated by the lower, dark gray band. The solid 
dark line, with gray uncertainty, and the light gray line are 
calculations using NLO fits to world g" data, evolved to our 
measured Q 2 range. 
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FIG. 2: Results for g% as a function of Q 2 are shown along 
with model calculations (see text). Data are shown with sta- 
tistical uncertainties only, with the systematic uncertainties 
indicated by the lower, dark gray band. The dark solid line, 
with gray uncertainty band, and the light gray line are calcu- 
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lations of g 2 using NLO fits to world (?™ data, evolved to Q 
our measured Q 2 . 



FIG. 3: Results for g% as a function of x for this experi- 
ment are shown along with previous measurements (see text). 
Data are shown with statistical uncertainties only. Results 
have been slightly displaced in x from the measured values 
to avoid overlap. The solid black curve shows gY W , without 
uncertainties, at Q 2 = 1.0 (GeV/c) 2 . 



Figure [5] shows our data for g% as a function of Q 2 
along with calculations ofg^ from the BB and AAC03e 
NLO analyses. The data are more than 5a above zero, 
and at lower Q 2 , show a systematic positive deviation 
from W i which we interpret as evidence for non-zero 
higher-twist contributions. In the OPE, the twist-2 and 
twist-3 contributions to gi enter at the same order in Q 2 , 
with additional higher-twist contributions suppressed by 
powers of 1/Q. Assuming these additional higher- twist 
contributions are small, we expect the quantity gi — g^ W 
to be constant as a function of Q 2 . When compared 
to g£ w from BB, a fit to our data gives 172 — 9i W = 
0.0262 ± 0.0043 ± 0.0080 ± 0.0099, with a reduced x 2 of 
1.4. The first two uncertainties are from the experien- 
tal statistical and systematic uncertainties and the third 
is from the uncertainty in the BB fit. Attempts to fit 
the data with functions that contain additional higher- 
twist dependence and/or logarithmic perturbative QCD 
corrections did not improve the quality of the fits. 

Also shown in Figure Elar e chiral soliton model calcula- 
tions from Weigel et al. |3lll32l | and Wakamatsu j3^|, and 
a bag model calculation of the higher-twist contribution 
from Stratmann |34|. combined with g^ w from BB. Our 
data indicate a positive contribution from higher-twist 
effects while the model calculations generally predict a 
negative contribution. 

Figure [3] shows the improved quality of our g% data 
(solid points) plotted versus x, as compared to older 
data (open squares) from the Stanford Linear Accelerator 
Center 0. Also shown are recent data (open triangles) 
from Jefferson Lab at higher x .23j. The data shown 
cover a range from Q 2 = 0.7 (GeV/c) 2 at the lowest x to 
20 (GeV/c) 2 at the highest x. For reference, gf w 
from BB is also included at constant Q 2 = 1.0 (GeV/c) 2 . 
In summary, we have made the first measurement of 
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the Q 2 -dependence of the gi spin structure function for 
the neutron at x ~ 0.2 in the range Q 2 = 0.57 — 1.34 
(GeV/c) 2 . With a factor of > 10 improvement in sta- 
tistical precision over previous measurements, these data 
allow for an accurate determination of the higher-twist 
contributions by direct comparison with the twist-2 g^ w 
prediction. Our results show a positive deviation from 
g^ w at lower-Q 2 , indicating a non-zero higher- twist con- 
tribution, and are inconsistent with model calculations, 
which generally predict a negative higher-twist contribu- 
tion. Precision data for gi obtained in this kinematic 
range showed good agreement with the NLO analyses of 
world data, indicating no significant higher-twist contri- 
butions within the uncertainties. 
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